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Abstract 

The exact Lyapunov dimension formula for the Lorenz system has been ana¬ 
lytically obtained first due to G.A.Leonov in 2002 under certain restrictions 
on parameters, permitting classical values. He used the construction tech¬ 
nique of special Lyapunov-type functions developed by him in 1991 year. 
Later it was shown that the consideration of larger class of Lyapunov-type 
functions permits proving the validity of this formula for all parameters of 
the system such that all the equilibria of the system are hyperbolically unsta¬ 
ble. In the present work it is proved the validity of the formula for Lyapunov 
dimension for a wider variety of parameters values, which include all param¬ 
eters satisfying the classical physical limitations. One of the motivation of 
this work is the possibility of computing a chaotic attractor in the Lorenz 
system in the case of one unstable and two stable equilibria. 

Keywords: Lorenz system, self-excited Lorenz attractor, Kaplan-Yorke 
dimension, Lyapunov dimension, Lyapunov exponents. 


1. Introduction 

The exact Lyapunov dimension formula for the Lorenz system has been 
analytically obtained first due to G.A.Leonov in 2002 [1] under certain re¬ 
strictions on parameters, permitting classical values. In his work it was used 
the technique of special Lyapunov-type functions, which had been created in 
1991 year [2] and then was developed in [3, 4]. Later in the works [5, 6, 7] 
it was shown that the consideration of a wider class of Lyapunov-type func- 
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tions allows to provide the validity of the formula for such parameters of the 
Lorenz system that all its equilibria are hyperbolically unstable. 

In this study it is proved the validity of the formula under classical restric¬ 
tions on the parameters. The motivation for this investigation is a numerical 
localization of chaotic attractor in the Lorenz system in the case of one un¬ 
stable and two stable equilibria [8, 9]. 

2. The Lorenz system 

Consider the classical Lorenz system suggested in the original work of 
Edward Lorenz [10]: 

{ X = a{y — x) 

y = rx — y — xz ( 1 ) 

z = —bz + xy. 

E. Lorenz obtained his system as a truncated model of thermal convection 
in a fluid layer. The parameters of this system are positive: 

a > 0, p > 0, 6 > 0, 

because of their physical meaning (e.g., b = 4(1 -|- a^)~^ is positive and 
bounded). 

Active study of the Lorenz system gave rise to the appearance and subse¬ 
quent consideration of various Lorenz-like systems (see, e.g., [11, 12, 13, 14]). 
A recent discussion of the equivalence of some Lorenz-like systems and the 
possibility of universal consideration of their behavior can be found, e.g. in 
[15, 16]. 

Since the system is dissipative and generates a dynamical system for Vt > 
0 (to verify this, it suffices to consider the Lyapunov function V{x, y, z) = 
|(a;^ + y"^ + {z — r — n)^); see, e.g., [10, 4]), it possesses a global attractor (a 
bounded closed invariant set, which is globally attractive) [17, 4]. 

For the Lorenz system, the following classical scenario of transition to 
chaos is known [8]. Suppose that a and b are hxed (we use the classical 
parameters a = 10, b = 8/3) and r varies. Then, as r increases, the phase 
space of the Lorenz system is subject to the following sequence of bifurcations. 
For 0 < r < 1, there is globally asymptotically stable zero equilibrium Sq. 
For r > 1, equilibrium Sq is a saddle and a pair of symmetric equilibria 
S'! 2 appears. For 1 < r < ^ 13.9, the separatrices ri 2 of equilibria 

So are attracted to the equilibria S'i,2- For r = Vh ^ 13.9, the separatrices 
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Fi 2 form two homoclinic trajectories of equilibria Sq (homoclinic butterfly). 
For Vh < r < Tc ^ 24.06, the separatrices Fi and F 2 tend to S 2 and Si, 
respectively. For rc < r < ra ^ 24.74, the equilibria S' 1,2 are stable and 
the separatrices F 12 may be attracted to a local chaotic attractor (see, e.g., 
[8, 9]). This attractor is self-excited^ and can be found using the standard 
computational procedure, i.e. by constructing a solution using initial data 
from a small neighborhood of zero equilibrium, observing how it is attracted, 
and visualizing the attractor (see Fig 1). 



Figure 1: (a) Numerical visualization of a self-excited Lorenz attrac¬ 

tor by using a trajectory with an initial point taken in the vicinity 
of the equilibrium So- (b), (c) Trajectories with the initial conditions 
(=F16.2899, =f 0.0601, 42.1214) tend to equilibria S' 2 , 1 . Parameters: r = 24.5, 
CT = 10, 6 = 8/3. 


For r > Ta, the equilibria S' 1^2 become unstable. The value r = 28 corre- 


^An oscillation can generally be easily numerically localized if the initial data from its 
open neighborhood in the phase space (with the exception of a minor set of points) lead 
to a long-term behavior that approaches the oscillation. Therefore, from a computational 
perspective, it is natural to suggest the following classification of attractors [18, 19, 20, 21], 
which is based on the simplicity of finding their basins of attraction in the phase space: An 
attractor is called a self-excited attractor if its basin of attraction intersects with any open 
neighborhood of an equilibrium, otherwise it is called a hidden attractor [18, 19, 20, 21]. 
Up to now in such Lorenz-like systems as Lorenz, Chen, Lu and Tigan systems only 
self-excited chaotic attractors were found. In such Lorenz-like systems as Glukhovsky- 
Dolghansky and Rabinovich systems both self-excited and hidden attractors can be found 
[22, 23, 24]. Recent examples of hidden attractors can be found in The European Physical 
Journal Special Topics: Multistability: Uncovering Hidden Attractors, 2015 (see [25, 26, 
27, 28, 29, 30, 31, 32, 33, 34, 35, 36]). Note that while coexisting self-excited attractors 
can be found by the standard computational procedure, there is no regular way to predict 
the existence or coexistence of hidden attractors. 
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spends to the classical self-excited local attractor (see Fig. 2). 



(a) Initial data near the equilib- (b) Initial data near the equilib- (c) Initial data near the equilib¬ 
rium So rium Si rium S2 


Figure 2: Numerical visualization of the classical self-excited local attractor 
in the Lorenz system by using the trajectories that start in small neighbor¬ 
hoods of the unstable equilibria S'o,i, 2 - Here the separation of the trajectory 
into transition process (green) and approximation of attractor (blue) is rough. 


3. Lyapunov dimension of attractors 

Consider a dynamical system 


x = f(x), (2) 

where f : —)■ M” is a smooth vector-function. Since the initial time is not 
important for dynamical systems, without loss of generality, we consider a 
solution x(t,xo) : a:(0,a:o) = xq. The linearized system along the solution 
x(t, xq) is as follows 


u = J(x(t, Xq)) u, u G M", t G M. 


( 3 ) 


where 

dfijx) : 
dxj l^=^(bxo) 

is the (nxn) Jacobian matrix evaluated along the trajectory x(t, xq) of system 
(2). A fundamental matrix X(t, xq) of linearized system (3) is dehned by the 
variational equation 


J(x(t,Xo)) = 


X(t,xo) = J(x(f,xo))X(f,xo). (4) 
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We usually set X(0,xo) = In, where is the identity matrix. Then u(t, Uq) = 
X(t, xo)uo. In the general case, u(t, uq) = X(t, xo)X“^(0, xo)uo. Note that if 
a solution of nonlinear system (2) is known, then we have 


X{t,xo) 


<9x(t,xo) 

<9xo 


Consider the transformation of the unit ball B into the ellipsoid X{t, xo)B 
and the exponential growth rates of its principal semiaxes lengths. The 
principal semiaxes of the ellipsoid X{t, Xq)B coincides with singular values of 
the matrix X(t, Xq), which are defined as the square roots of the eigenvalues of 
matrix X(t,Xo)*X(t,Xo). Let (Ti(X(t,xo)) > ■■■ > (T„(X(t,xo)) > 0 denote 
the singular values of the fundamental matrix X(t,xo) reordered for each t. 
Introduce the operator T’(-) = Iln| ■ |, where | ■ | is the Euclidian norm. 
Dehne the decreasing sequence (for all considered t) of Lyapunov exponents 
functions LEj(t,a;o) = X{ai(t,Xo)) 


LEi(t,Xo) > LE 2 (t,Xo) > ... > LE„(t,a;o). 


Definition 1. [37] The Lyapunov exponents (LEs) at the point xq are the 
numbers (or the symbols ±oo) defined as 

LEj(xo) = limsupLEj(t, xo) = limsup - ln(Tj(X(t, xq)). (5) 

t—>00 t—KX) t 

LEs are commonly used^ in the theory of dynamical systems and dimen¬ 
sion theory [46, 47, 48, 49, 50, 4, 51]. 


^Two widely used definitions of Lyapunov exponents are the upper bounds of the 
exponential growth rate of the norms of linearized system solutions (LCEs) [38] and 
the upper bounds of the exponential growth rate of the singular values of fundamen¬ 
tal matrix of linearized system (LEs) [37]. The LCEs [38] and LEs [37] are “often” 
equal, e.g. for a “typical” system that satisfies the conditions of the Oseledec theorem 
[37]. However, there are no effective rigorous analytical methods for checking the Os¬ 
eledec conditions for a given system [39, p.ll8] (a numerical approach is discussed in 
[40]). For particular system, LCEs and LEs may be different. For example, for the 

1 

0 1 

LCEi = max(limsupT[g(t)],limsupT[(/“^(t)]),LCE 2 = 0; 

>- + 00 

LEi ,2 = max,min(limsupT[(/(t)],limsupT[ 5 “^(t)]). Note also that positive largest 

t —^~|”00 t —^“|“00 

LCE or LE, computed via the linearization of the system along a trajectory, does not 


fundamental matrix X(t) = 


we have the following ordered values: 
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Remark 1. The LEs are independent of the choice of fundamental matrix 
at the point xq [45, 52]. 

Consider the largest integer j(t,Xo) G {1, ..,n} such that 

LE°(t, xo) + ... + LE](^^^^^^(t,xo) > 0, 


LE"(t,a;o) + ... + Xq) 

I j{t,xo)+li^^ ^o) I 


< 1 . 


Following [53, 49], introduce the following dehnition: the function LD(t, xq ) 
0 if LE°(f,a;o) < 0 and LD(f, xq) = n if Y[h=i LE°(C2^o) > 0, otherwise 


LD(f,a;o) = j{t,XQ) + 


LE]'(f,a;o) + ... + Xq) 


LE' 


'j{t,xo)+l{^i ^o) I 


( 6 ) 


Definition 2. A local Lyapunov dimension at the point xq is as follows 

LD(a;o) = linisupLD(f, xq). 

t —^~1“CXD 

The Lyapunov dimension of invariant compact set K of dynamical system 
is dehned by the relation 


dim^iL = sup LD(a;o) = sup limsupLD(t, xq) 

xq'^K xq'^K t^+oo 


Note that, from an applications perspective, an important property of 
the Lyapunov dimension is the chain of inequalities [49, 54, 4] 

dimT’iL ^ dim/f iL ^ dimjT’iL ^ dim^ K, (7) 

where dim^, dimj^ iL, and dim^^iL are topological, Hausdorff, and fractal 
dimensions of K, respectively. 


necessary imply instability or chaos, because for non-regular linearization there are well- 
known Perron effects of Lyapunov exponent sign reversal [41, 42, 43]. Therefore, in general, 
for the computation of the Lyapunov dimension of attractor we have to consider a grid of 
points on the attractor and corresponding local Lyapunov dimensions [44]. More detailed 
discussion and examples can be found in [45, 41]. 
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4. Estimation of Lyapunov dimension by Lyapunov functions 

Along with commonly used numerical methods for estimating and com¬ 
puting the Lyapunov dimension (see, e.g., [44, 24]), there is an analytical 
approach that was proposed by G.A.Leonov [2, 3, 4, 55, 5, 15]. It is based 
on the direct Lyapunov method and uses Lyapunov-like functions. The ad¬ 
vantage of this method is that it allows to estimate the Lyapunov dimension 
of invariant set without numerical localization of the set. This is especially 
important for the systems with hidden attractors when numerical hnding of 
all local attractors may be a challenging task [24]. 

Since LEs and LD are invariant under the linear changes of variables (see, 
e.g., [45]), we can apply the linear variable change y = Sx with a nonsingular 
n X n-matrix S. Then system (2) is transformed into the system 

Y = s± = s'f(S'"V) = f(y)- 

Consider the linearization along the corresponding solution y(t, yo) 
that is, 

V = J(y(t,yo)) V, V e 
Here the Jacobian matrix is as follows 

= S J{x{t,xo))S-^ (9) 

and the corresponding fundamental matrix satishes Y(t,yo) = SX(t,XQ). 

For simplicity, let J(x) = J(x(t,xo)). Suppose that 

Ai(x, S') ^ ■■■ ^ An(x, S') are eigenvalues of the symmetrized Jacobian 
matrix (9) 

l(SJ(x)S-' + (SJ(x)S-')-). (10) 

Theorem 1 ([5]). Given an integer j G [l,n] and s G [0,1], suppose that 
there are a continuously differentiable scalar function -d : M” —)■ M and a 
nonsingular matrix S such that 

Ai(x, S') J-■ ■ ■ J-Aj(x, S') J-sAj+i(x, S') J--d(x) < 0, VxgJL. (11) 
Then dimj) K ^ j + s. 

Here d is the derivative of d with respect to the vector held f: 

d(x) = (grad(d))T(x). 


= Sx{t,S %), 
( 8 ) 
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Theorem 2 ([2, 3, 4, 5]). Assume that there are a continuously differentiable 
scalar function A and a nonsingular matrix S such that 

Ai(x, S) + A 2 (x, S) + ^9(x) <0, Vx e ( 12 ) 

Then any solution of system (2) bounded on [0, +cx)) tends to a certain equi¬ 
librium as t ^ +CX). 


5. Main result: Lyapunov dimension of the global Lorenz attractor 

By Theorems 1 and 2, for the Lorenz system we can obtain the following 
result. 


Theorem 3. Assume that the following inequalities 

r-l>0, (13) 

+ a — 1)^ — Aa{h -\- ah — b‘^) 

3(j2 

are satisfied. Let one of the following two conditions be satisfied: 

a. 




a'^{r — 1)(6 — 4) < Aa{ah + h — h"^) — h{h + a — 1)^; (15) 

b. there are two distinct real roots of equations 

(2a — & + 7 )^ (b(b + a — 1)^ — 4a(ab + 6 — 6 ^) + a^(r — l)(b — 4)) + 

+ 4 & 7 ((T + 1) (6(5 + a — 1)^ — 4a(ab 4- b — b^) — 3a^(r — 1)) = 0 

(16) 


and 

a‘^(r — l)(b — 4) > 4a(ab + b — b"^) — b(b 4- a — 1)^ 

7 ("") > 0 

where is a greater root of equation (16). 

In this case 
1 . if 

(b — a)(b — 1) < ar < (b 4- l)(b + a), (18) 

then any bounded on [ 0 ; +cx)) solution of system ( 1 ) tends to a certain 
equilibrium as t ^ +cx). 




ar > {b+ l){b + a), 


(19) 


2. If 


then 


diiiii K < 3 — 


2(cr + 6+1) 

a + 1 + a/(ct^-TP~+4^’ 


where K is a bounded invariant set. 


( 20 ) 


For numerical experiments (see, e.g., [56, 57]) it is known that the Lya¬ 
punov dimension of any invariant compact set of the Lorenz system is bounded 
from above by the local Lyapunov dimension of the zero equilibrium. 


Lemma 1. If for the parameters of system (1) inequality (19) is valid, then 


dim2,(0, 0, 0) 


2((j + 6+1) 

a + 1 + A/(cr — 1)^ + Ara 


Comparing estimation (20) with the local dimension in the origin, we 
obtain a dimension formula for a global Lorenz attractor. 


Theorem 4. Let K be a global attraetor of classical Lorenz system (1) and 
conditions (13)-(17) of Theorem 3 he hold true. 

In this case if 

or > (6 + 1)(6 + a), 


then 


dim^ K 


2((T + 6 + 1) 

a + 1 + a/((T — 1)^ + dor 


( 21 ) 


Remark 2. It can he easily checked numerically that if all three equilibria 
are hyperbolic, then the conditions of Theorem 4 (see Fig. 3b) are satisfied. 
For example, for the standard parameters cr = 10 and 6 = | the formula (21) 
is valid for r > ^. 


5.1. Numerical analysis 

We perform a numerical analysis of parameter domain, which does not 
satisfy the conditions of Theorem 4. 

1. Recall that for 0 < r < 1 the solutions of system tend to a unique 
equilibria (see., e.g., [4]). 

2. Consider r > 1. In the case of conditions (13)-(19) the following 
assertion is valid. 
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(a) r = 10 



Figure 3: Numerical simulation of domains for the fixed values of parameter 
r = 10,28,100 


Lemma 2. For all r > 1, a > 7, and 0 < 6 < 4 the conditions of Theorem 
3 are satisfied. 

Remark 3. For 0<6<4 — e, ct>7 and r > r^, where r* = r*(e), the 
conditions of Theorem 3 are satisfied. 

Thus, for each fixed r > 1 it remains to check numerically the behavior 
of system with the parameters from the bounded domain [0 < 6 < 4, 0 < 
a < 7]. For the considered increasing sequence rj. > 1, numerical simulation 
shows the lack of chaotic attractor (the trajectories with initial data from 
the corresponding compact absorbing set [58] have been simulated). 

Appendix 

Proof. (Theorem 1) 

For the proof we make use of Theorems 1, 2. Consider a solution of 
system (1) x= {x{t),y{t), z(t)) = {x,y,z). 

I. In the case of classical Lorenz system the matrix J has the following 
form 



Following [4], for the condition 

or + (a — 6) (6 — 1) >0 (22) 


10 











we introduce a matrix 


/-p-i 0 0\ 

1 Oh 

Vo 0 1/ 


where 


Then 


^ ^/oTTJa^^b){b^^ 

j h-a-l 0\ 

--p+pz -b - x \ 
\-p {y + ^x) X -b) 


( 23 ) 


We hnd the eigenvalues of the matrix ^{SJS ^ + {SJS ^)*). The character¬ 
istic polynomial of the matrix takes the form: 


(A+ 6) 


A^ -h (a -h 1)A b{a 


6 + 1 ) 




1) p 

—x+^y 


2 


This implies that the eigenvalues A* = \ i { x , y , z , S),i = 1,2,3 of the 
matrix ^{SJS~^ + {SJS~^)*) are the numbers: 

A2 = -b, 

(a-2b+lf+(^j-pzy+ f/(y + ^xj" \ 

Such a choice of the matrix S provides a simple form of eigenvalues. 


-^1,3 — 


a 


±t±i 
2 2 


II. We shall show that Ai > A 2 > \ 3 ,\/x,y,z: 

-{a + 1) + {a - 2b + 1)^ + (^^ - pz'^ + 


> -b> 


> 


-{a + 1) - ^{a - 2b + 1)^ + (^^ - pz'^ + p"^ {y + 


2 
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> a — 2b + 1 > 


> - 

a — 2b + 1\ 

At the left there is a sum of two total squares, i.e. the inequality is always 
satished. 

Thus, Ai > A 2 > \ 3 ,\/x,y,z. 

III. To apply Theorem 1 and Theorem 2, for s G [0,1), we consider and 
transform the relation 2(Ai + A 2 + sAs): 

2(Ai + A 2 + 3 X 3 ) = —((T + 26 + 1) — s((j + 1)+ 

+ (1 - s) (a - 26 + 1)^ + {y + 

— ((T + 26 + 1) — s((T + 1) + 

+ (1 - s) (a - 26 + 1)^ + ^ - Aaz + p^z"^ + (y + -—-x 

= —((T + 26+1) — s {(7 + 1) + 

{a — 1)^ + Aar — Aaz + p‘^z‘^ + p^ 

(24) 

We can use the inequality 

Vk + I < Vk H-^ (25) 

~ 2Vk 



+ il-s) 
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for /c + / > 0, /c > 0. Consider 


h — 1 

k + I = {a — ly + Aar — Aaz + [y -\ - x 


and 


k = {a — ly + Aar. 

According to (24) the following equation 


k + I = {a — ly + Aar — Aaz + p'^z^ + p ( y H- x ) = 


h-l 


2a 


= (a — 2& + 1)"^ +- pz\ + p^ y -\ - X > 0 


P 


b-1 


a 


is valid. By the condition of Theorem 3 we have 

ra + {a — b){b — 1) > 0 ^ ra > —ab + a + ^ — b. 

Then 

k = {a — ly + Aar > {a — 1)^ — Aab + Aa + Ab"^ — Ab = {a — 2b + 1)^ > 0. 

Applying the inequality (25) with k = (a — 1)^ + Aar and / = —Aaz + 
p^z"^ + p^ {y + ^x) to the relation (24), we obtain 

2(Ai + A 2 + <sA 3 ) < —[a + 25 + 1) — s((T + !) + (! — 'S)[(cr — 1)^ + Aary-P 


+ 


2(1-5) 


{a — ly + Aar] 2 


p^z‘^ p^ 


A A 


b-1 

az+'— + ^iy+ -—X 


a 


Introduce the function 


0{x,y,z) = 


{1 - s)V{x,y,z) 
[(a — ly + 4(Tr]5 


(26) 


(27) 


where 


V{x,y,z) = ^ 4 X^ + {-a'yi+^ 3 )y‘^+j 3 z‘^ + ^jix‘^-^ix^z-jij 2 xy-^z. (28) 
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We choose the running parameters 7 i, 72 , 73,74 in such a way that 




h-l 


i? ;=-CT^ + ^—+ — + ——xj +V <0 yx,y,z, (29) 

where V{x) is a derivative with respect to system (1). Then, using (26), we 
obtain 

2(Ai + A 2 + <sA 3 )+ 20 < —((T + 26 + 1) — s(cr + l) + (l — >s)[(cr — l)^ + 4(Tr]2, (30) 

i.e. we estimated 2(Ai + A 2 + sX^) + 29 by the relation, which depends on 
the parameters of system ( 1 ) and is independent of x,y,z. 

IV. We perform the analysis of R, choosing the rnnning parameters 
11 , 12 ,Is, 14 : in snch a way that (29) is valid. Taking into acconnt (28), we 
obtain 

r, , I,., , f “ 1 )^ ^ ...2 , 


R = -'fix + (7172 + 2 (T 7 i + kyi)x z + 


- ''"7172 - 2(774 X + 


pHb-1) 


+ I 1 I 2 - 2 ra 7 i + 2773 + 20-74 + 0-7172 - - 1 xy+ 


+ 20-71 - 273 - 0-7172 + ^ ^ - 2073 = 


= Aix‘^ + A2 x‘^z + A'iZ^ + Bix^ + B2xy + B^^y"^, 

where 

P^ 

Ai = -71, A2 = 71(72 + 2o- + 6), A3 = — - 2673 
p^(b — 1 )^ 

Bi = - ^ -r 7 i 72 - 20-74 = Cl - 20-74 

B2 = —— + I1I2 - 2ro-7i + 2r73 + 0-7172 - ^ + 2o'74 = C2 + 20-74 

2 o- b 

p^ 

Bs = 20-71 - 273 - 0-7172 + —. 


IV. a. Under the conditions Ai 7 ^ 0 and B^ ^ Owe transform the relation 
for R: 


R = Ai{x^ + —z] +BAy+—x] + 


451^3 -Bl 2 , 4 A 1 A 3 - A 2 
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Then 


Ai<0 

B3 < 0 

4^1^3 ~ ^2 — 0 
?2 


> i? < 0 Vx, y, z 


ABiB^ — B2 ^ 0 J 

Now we analyze the signs of four addends in (31) 

( 1 ) 

Ai = —71 < 0 77 7 i > 0 

( 2 ) 


( 3 ) 


P 

-63 < 0 77 273 > 2(T7i - (77172 + — 


0 < AAiAs - ^2 = -471 ( ^ - 2673 ) - 71^(72 + 2(7 + 6 ) 


Since according to (33), 71 > 0, from (35) it follows that 
- 71(72 + 2(7 + 6 )^ - + 8673 > 0 


(32) 


(33) 

(34) 


(35) 


77 

4-01-03 — -02 = 4 - 03 (C'i — 2 ( 774 ) — {C 2 + 2 ( 774 )^ = 
= 4-03Ci — SaBs'j^ — cl — 4(7(0274 — 4(7^74 = 

= + 474a(-253 - C 2 ) + 453^1 - Cl 


(36) 


The relation 4 - 01-03 — -0| is a quadratic polynomial in 74 with a negative 
coefficient of 7 |. Therefore for the inequality 4 - 01-03 — -0| > 0 to be satisfied, 
it is necessary that the corresponding quadratic equation has a real root, i.e. 
a positive discriminant: 


= IQa^ABl + Cl + 4 - 03 ^ 2 ) + IQa^AB^Ci - Cl) = 

= lQCla‘^+Ua‘^Bl+Ua‘^B3C2+Ua‘^B3Ci-lQa‘^Cl = 64a^-03(-03+C'2+C'i). 
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Since according to (32) the condition ^3 < 0 is satisfied, we obtain > 
0 -v^ -83 + Cl + C 2 < 0. We have 


-83 + Cl + C 2 


2 (T 7 i - 273 - (J7i72 + 


C 

4 


+ 


p‘^{b — 1 )^ 
4(7^ 


— 1) a 

-^7172 H- 7 -h 7 i 72 - 2ra7i + 2773 + 0-7172 - t = 

2 o- b 

p^ib — 1)^ 

= 73(2r - 2) - 7 i(-2o- + 2ro-) - 7172(0- + r- o--l) + ^ + ^ 4^2 + 


p ^(6 — 1 ) o- 


P 


2^ - - = 273(7 - 1) - 2^71 (r - 1) - 7 i72(?’ - 1) + ^ + 

4o-2 2 o- b 


Thus, for 


273(7 - 1) < 20-71(7 - 1) + 7 i 72 (?’ - 1) - 

p‘^{b — 1 )^ p‘^{b — 1 ) a 

4o-^ 2o- 6 


(37) 


there exists 74 such that ABiB^ — > 0 . 

Note that a part of the right-hand side of the inequality in (37) can be 
transform in the following way 


P^ , P\h-lf p\b-l) 

4 4o-2 2o- 

_ ^ {b-lf cT{b-l) \ _ p\a + b-lf 

o-^\4~^ 4 2 ) 4o-^ 


(38) 


Taking into account (33), (34), (36), (37), (38), conditions (32) become 


7i > 0 

P^ 

273 > 20-71 - 0-7172 -h — 

2(7 - 1)73 < 20 - 71(7 - 1 ) - 1 - 7i72(»^ - 1 ) 


> ^ -R < 0 Vx, y, z. 


p\b + a-ir a 
4o-2 b ^ 

( 39 ) 
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For obtaining ( 31 ) it is assumed that Ai ^ 0 ,B^ ^ 0 . Let us analyze the 
cases Ai = 0 and = 0 . 


IV. b. Consider B^ = 0 . In this case we have 

R = Aix"^ + A2 x‘^z + Asz^ + Bix^ + B2xy, 

i.e. it is impossible to choose 71,72,73,74 in such a way that i? < 0 is valid 
'ix,y,z. 


IV. c. Consider the case 71 = Ai = 0 and -B3 ^ 0 . Then 

A2 = 71(72 + 2(7 + 6 ) = 0 


and 


-R — + Six + B2xy + — 

2 , n / , B2 ^ 4 :BiBs 


— A3Z + Bs [ y + 


25, 


-X 


+ 


4:B:, 


Bl 2 

—X . 


In this case 

dls <0 1 

53 <0 >^-R <0 \/x,y,z. ( 40 ) 

45i53 - ^2 > 0 j 

The second and third conditions are similar to the second and fourth condi¬ 
tions in ( 32 ). Consequently it remains to consider 

2 2 

0 > ^3 = ^ - 2673 77 273 > 

4 46 

The latter inequality, obtained under the assumption Ai = 0 , coincides with 
condition ( 34 ) if in ( 34 ) we consider 71 = 0 . Thus, the conditions on the 
running parameters 71,72,73,74, obtained under the condition Ai = 0, ^3 7^ 
0, can be joined with those, obtained under the condition Ai 7^ 0, ^3 7^ 0, 

i.e. 


7i > 0 


273 > 2(T7i - (77172 -h 


4 


2(r - 1)73 < 2(771 (r - 1) -1- 7i72(»^ - 1) 


> ^ i? < 0 Vx, y, z. 


p\b + a-ir a 
4(7^ 6 ^ 

( 41 ) 
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V. For r — 1 > 0 conditions (41) can be transformed in the following way: 


7i > 0, 

273 > ^ + 2(771 - (77 i 72, 

P =(6 + <t- 1 )" 

273 < 2cr7i + 7 i 72-r^TT- 73 -1- 


(42) 


(7 


4(7^(r — 1) b{r — 1) 
For the existence of 73 it is necessary and sufficient that 


+ (7 — 1)^ 

— + 7 i (2(7 - (772) < 7 i (2(7 + 72) - 


+ 


(7 


4(7^(r — 1) b{r — 1) ’ 
+ a - If 


a 


^ + 1(2,3 + (, + 73)^ < 71(2.^ + 7.) - _ 1)' + ■ 

We now analyze the obtained inequalities. 

V. a. Consider the first inequality from system (43) for different 72. 
Let be 72 = 0. Then the inequality is equivalent to 


(43) 


p^{b + a — lf (7 
4 Aa‘^{r — 1) bfr — 1) 

If 72 > 0, then we obtain the condition 
I / ^2 ^2 


< 0 . 


7i > 


l^p\b + a-iy 


a 


72((7 + 1) \4 4a‘^{r — 1) b{r — 1) 


and if 72 < 0, then 


7i < 


p^ ^ p^{b + a — if 


a 


72((7 + 1) \4 4(72(r — 1) 6(r — 1) 

In the latter case it is required that 

-. \ O 

a 


p^ p^fb + (7 — 1 )^ 

4 4(7^(r — 1) b{r — 1) 


< 0 


(44) 
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since according to (41) we have 71 > 0. Thus, since the cases of 72 = 0 and 
72 < 0 impose the same condition on the parameters of the system, they can 
be joined. 


V. b. From the second inequality of system (43) we obtain condition for 


7i: 


(2a-b + 72 )^ 
46 


p2(b + a-lf 


a 


46 4(T2(r — 1 ) b(r — 1) 


If 72 = 6 — 2a, then it is required that 

p^(h + a — ly a 
46 Aa‘^(r — 1) b(r — 1) 


< 0. 


If 72 7 ^ 6 — 2a, then we obtain the following condition for 71 

46 


7i < 


^ p^(b + a — 1 )^ 


a 


(2(T — 6 + 72 )^ \46 4a‘^(r — 1) 6 (r — 1) 

According to (41), it is valid 71 > 0. Therefore it is required the following 


p^ p^(b + a — ly a 
46 Aa‘^(r — 1) 6 (r — 1) 


< 0. 


(45) 


Thus, inequality ( 45 ) holds true also for 72 = 6 — 2 (T and for 72 7^ 6 — 2a. 

Condition (45) must be satished for any sign of 72 . If, in addition, con¬ 
dition (44) is also satished, then in function (28) we can take 72 < 0 and in 
this case there exist 71 , 73,74 such that the relation 

[ 6 - 1 \ ^ 

R=—az+^^ — ^”4 - ^x,y,z 


is valid. If inequality (44) is not satished, then in function (28) it must 
be considered 72 > 0 and, except for condition (45), it is necessary to hnd 
conditions for the existence of 71 . 

In this case we obtain the following family: 
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' (p^ , Ab + ^-if _ 

I 4b 4(T^(r — 1 ) b{r — 1 ) “ ’ 

1 , p\b + a-ir _^ 

I 4 4a'^{r — 1) b{r — 1) 

r P^ , p\h + a-lf _ 

4b 4a‘^{r — 1) b{r — 1) “ ’ 

72 > 0,72 7 ^ & - 2 ( 7 , 


< 1 (p^ p\b + <j-iy a \ 

72(0- + 1) \ 4 4a‘^{r — 1) b{r — 1)/ ^ 

4b f p^{b + a — lY a 

^ (2(7 — 6 + 72)^ \46 4(7^(r — 1) b(r — 1) 


( 46 ) 


V. c. 

Consider the first condition of family (46): 

p^ p^{b + a-lY cr 
4b 4(7^(r — 1) b(r — 1) 

p^ p^{b + a — lY O' 

4 4a‘^{r — 1) b{r — 1) 


< 0 , 


< 0 . 


Substituting the value p from (23), we obtain 

(h + a-lf 


+ 


46((r — 1)(7 — b‘^ + b + ab) 4(r — l)((r — l)a — b‘^ + b + ab) 


(7" 


+ 


(6 + (7 - 1)" 


4((r — 1)(7 — b'^ + b + ob) 4{r — l)((r — l)a — b'^ + b + ob) 


■77 


r — 1 > 


6(6 + (7 — 1)^ — 4(7(6 + (76 — 6^) 
3(7^ 


(7^(r — 1)(6 — 4) < 4a{ab + 6 — 6^) — 6(6 + (7 — 1)^ 


(7 

b{r — 1) 

(7 

6(r — 1) 


< 0 , 

< 0 . 


20 



V. d. Consider the second condition of family (46): 


' , p\h + a-lY _ 

46 4(T2(r — 1) b{r — 1) “ ’ 

72 > 0 , 72 7 ^ 6 - 2a, 

< 1 p^(6 + a-l)^ _ g \ (47) 

72(g + 1) \ 4 4g2(r — 1) b{r — 1)/ 

46 P p^(6 + g — 1 )^ a \ 

^ (2g — 6 + 72)^ V46 4g2(r — 1) 6(r — 1) / ’ 

From Item V. c. it is known that the first inequality of system is equiv¬ 
alent to the relation 


r — 1 > 


6(6 -I- g — 1)^ — 4g(6 -|- g6 — 6^) 


Next we transform the inequality 


1 (I , p\b + a-ir _ a \ 
72(g -I- 1) \ 4 4g2(r — 1) 6(r — 1)/ 

46 / p^(6-|-g —1)^ g 

(2g — b + 72)^ \46 4g2(r — 1) b{r — 1) 


\2 f P^ p‘^{b + a — iy a \ 

( 2 <t -6 + 72 ) (7 + ^^5^737-- iyvT)) + 

+4M- + 1)(d + dh±^_ - )<o 

\46 4g2(r — 1) b{r — 1)/ 


(2g - 6 -h 72 )^ [ 6 ( 6 -h g - 1)^ - 4g(g6-h 6 - 6 ^)-h g^(r - 1)(6 - 4)] 

-|- 4672 (g -I- 1) ( 6(6 -I- g — 1)^ — 4g(g6 -|- 6 — 6 ^) — 3g^(r — 1)) < 0. 

The left-hand side of the inequality is a quadratic polynomial in 72. If 
coefficient of 7I is negative, i.e. 

6(6 -I- g — 1)^ — 4g(g6 -|- 6 — 6 ^) -|- g^(r — 1)(6 — 4) < 0, 


then there exists 72 > 0 satisfying inequality (48). 
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If the coefficient is equal to 0, i.e. 

b{b + a — 1)^ — Aa{ab + b — b'^) + a^{r — 1)(6 — 4) = 0, 
then the relation 

4672(0- + 1) [6(6 + O' — 1)^ — 4o-(o-6 + 6 — 6^) — 3o-^(r — 1)] <0 

must be valid. By (47), 72 > 0 is satisfied and 6(6 + o' — 1)^ — 4o-(6 + o-6 — 
6^) — 3o-^(r — 1) < 0, Then the inequality 

4672(0- + 1) [6(6 + 0- — 1)^ — 4o-(o-6 + 6 — 6^) — 3o-^(r — 1)] <0 

is satished. 

If the coefficient is positive, i.e. 

6(6 + 0- — 1)^ — 4o-(o-6 + 6 — 6^) + a^ir — 1)(6 — 4) > 0, 

then for the existence of 72 > 0, satisfying (48), it is necessary and sufficient 
that there exist two different real roots of the equation 

(2o-- 6 + 72)^ [6(6 + 0-- 1)^ - 4o-(o-6 + 6 - 6^) + o-^(r - 1)(6 - 4)] + 
+4672(0- + 1) (6(6 + 0- — 1)^ — 4o-(o-6 + 6 — 6^) — 3o-^(r — 1)) = 0 

and the largest root must be positive. This condition corresponds to the 
second condition, stated in the theorem. 

V. e. Thus, there are obtained conditions (14)-(17) for which 

f b-l V 

R = —az + ^ X \ ^ -^ ^ ^ 

Consequently, it is valid (30): 

2(Ai + A 2 + <sA 3)+20 < —( 0 - + 26 + 1) — s(o- + 1) + (1 — s) [( 0 - — 1)^ + 4o-r] 2 . 

VI. To apply Theorem 2, we consider inequality (30) for s = 0; 

2(Ai + A 2 + ^) ^ —( 0 - + 26 + 1) + [( 0 - — 1)^ + 4 o-r] 2 . 
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Then we find the conditions for the parameters of system a, b, r for which 
2(Ai + A 2 + 0) < 0. We have 

— ((T + 26 + 1) + [((T —l)^ + 4(Tr]^ < 0 [{a — + Aar]^ < {a + 2b + l). (50) 

The parameters of system a,b are positive numbers, i.e. <7 + 26+1 > 0. 
Therefore after squaring two sides of inequality (50) we obtain 

(7^ — 2(7 + 1 + 4(77 < (7^ + 1 + 46^ + 2(7 + 4(76 + 46 


4(77 < 46^ + 4(7 + 4(76 + 46 

(77 < (6 + 1)(6 + a ) 

Thus, taking into account inequality (22), we obtain condition (18) for which 

2(Ai + A 2 + ^) <0, 

and Theorem 2 can be applied, i.e. any bounded on [0; +cx)) solution of 
system ( 1 ) tends to a certain equilibrium as f —)■ +cx). 

VII. To apply Theorem 1, we consider s ^ 0 and (30): 

2(Ai + A 2 + sAs + 0) < -((7 + 26 + 1) - s(a + 1) + (1 - s )[((7 - 1)^ + Aar]^ = 

= — s((7 + 1 + [((7 — 1)^ + 4(77] 2 ) — (a + 26 + 1) + [((7 — 1)^ + 4(77] 2 . 

We find s such that the relation 2(Ai + A 2 + sX^ + 0) < 0 is valid: 

— s((7 + 1 + [((7 — 1)^ + 4 (77] 2 ) — ((7 + 26 + 1) + [((7 — 1)^ + 4 (77] 2 < 0 


— {a + 26 + 1) + [((7 — 1)^ + 4(77]2 < s((j + 1 + [((7 — 1)^ + 4(77]2). (51) 

The parameter <7 of system (1) is a positive number. Therefore the coef¬ 
ficient of s, equal to <7 + 1 + [(a — 1)^ + 4 ( 77 ] 2 , is greater than 0. In this case 
if inequality (51) is divided by a + 1 + [(a — 1)^ + 4 ( 77 ] 2 , then the sign is not 
changed, i.e. we have 


— ((7 + 26 + 1 ) + [((7 — 1 )^ + 4 ( 77 ] 2 
(7 + 1 + [((7 — 1)2 + Aar] 5 


(52) 
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Thus, in the case when (14)-(17) is satished and s satishes (52) we have 
2(Ai + A 2 + SA 3 + 0 ) < 0 . 

Under the hypothesis of Theorem 1 we have s G [0,1). However accord¬ 
ing to (52) the relation s > Sq must be valid. The case s = 0 is already 
considered. Obviously a -|- 1 -|- [(n — 1 )^ -|- dor] 2 > 0. Therefore for sq > 0 to 
be valid, it is required that 

— {a + 2b+ 1) + [(a — 1 )^ -|- darjs > 0 


[(a — 1)^ -I- 4(Tr]2 > a -I- 25 -I- 1. (53) 

The relation a -|- 26 -|- 1 > 0 is always satisfied since the parameters a, b of 
system (1) are positive numbers. Thus, after squaring inequality (53) we 
obtain 

— 2a + 1 + Aar > -|- 1 -|- 46^ + 2a + Aab + Ab 

or > (6 -I- 1 )(6 -I- a). (54) 

Next we hnd a condition for which Sq < 1: 

-{a+ 2b+l)+ [{a-ly+ Aar]^ 
a -|- 1 -|- [(a — 1)^ -|- Aar] 2 

— {a -|- 26 -|- 1) 4- [(a — 1)^ -|- dor] 2 < a -|- 1 -|- [(a — 1)^ -|- dor] 2 

0 < a -I-1 -I- 6 

The latter inequality is always valid. The parameters of system (1) are 
positive numbers and therefore 


(6 + l){b + a) > (6 - a )(6 - 1 ) 


is always satished, i.e. from conditions (22), (54) it remains only condition 
(54). Thus, it is obtained condition (19) under which the relation 2(Ai -|- 
A 2 + -sAa -|- 0) < 0 is satished and Theorem 1 can be applied. Consequently, 
dim^ K < 2 + s for s, satisfying (52) and (19). Thus, 


dimz, K <3 — 


2((T -|- 6 -|- 1) 

a -|- 1 -|- \/{a — ly + Aar 


□ 


24 



Proof. (Lemma 1) For positive parameters we have 

{b + l){b + a) > {b - a){b - 1). 

Consequently if condition (19) is satisfied, then the inequality 

ar > {b— l){b — a) 

is satished too. First we transform this inequality 

ra + {a — b){b — 1) > 0 6 ^ — b{a + 1 ) + la — ra < 0 

bi<b < 62 , (55) 

where 61,2 = 

When considered the linearized system along the zero solution the corre¬ 
sponding matrix J(x) = J is constant: 

/—a a 0 \ 

J = I r —1 0 ) • 

\ 0 0 - 6 / 

Its eigenvalues are real and have the form: 

Ai = -^(a -h 1 - A/(a - 1)2 + Ara) 

\2 = -b 

As = -^(a -h 1 -h A/(a - 1)2 -h 4ra). 

Then LEj(0, 0, 0) = Aj, i = 1, 2, 3. 

Condition (55) is valid: 

b 2 > b > bi (a - 1 )^ 4ra) > -b > (a - 1 )^ 4ra). 

Thus, LEi(0,0,0) > LE 2 ( 0 , 0 , 0 ) > LE 3 ( 0 , 0 , 0 ). 

We shall show that if (19) is satished and j(t, (0,0,0)) = 2, then the 
dehnition of ( 6 ) satishes all requirements of the dehnition of local Lyapunov 
dimension. 
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( 1 ) 

LEi(0, 0, 0) + LE2(0, 0,0) = + 1 - a/((J - 1)^ + Ara) - b = 

= -^(o- + 26 + 1) + ^a/((T - 1)2 + ira. 

The relation LEi(0, 0, 0) + LE2(0, 0, 0) > 0 must be valid: 

-^((T + 26 + 1) + ^A/(a - 1)2 + 4ra > 0 
a/((J — 1)2 + 4ra > (j + 26 + 1 

We obtained inequality (53), which, as is known, is equivalent to inequality 
(19). 

(2) 

LE3(0, 0, 0) = -^(a + 1 + A/(cT^^Tj2”+4m) < -6. 

Parameter 6 of system (1) is positive, consequently LE3(0,0,0) < 0 is satis- 
hed. 

(3) 

LEi(0,0, 0) + LE2(0, 0, 0) _ -|(a + 1 - a/(ct - 1)2 + 4ra) - 6 _ 
|LE3(0,0,0)| I - 1((T + 1 + A/(cr - 1)2 + dra)! 

_ -(a + 26 + 1) + A/(a - 1)2 + 4ra 
(T + 1 + \/ (cT — 1)2 + Ara 

The expression obtained coincides with sq of inequality (52). The relation 

LEi(0, 0, 0) + LE2(0, 0) 0) _ ^ ^ 

|LE3(0,0,0)| 

must be valid. According to the proof of Theorem 3, we have Sq < 1 
a + 1 + 6 > 0. This is true since the parameters of system a, 6 are positive 
numbers. 

Thus, we have 

dim,(0,0,0) = . a + 2t+l-y(a-l)^ + 4m_,- 2{a + b+l) 

(J + 1 + a/((T — 1)2 + 4ra a + 1 + a/(o' — 1)2 + 4ra 

□ 
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Proof. (Lemma 2) 

I. We consider sufficient condition for the theorem to be valid: 
b{b + a — 1)^ — 4a{b + ab — 6^) 

^ (56) 

(T^(r — 1)(6 — 4) < Aa{ab + b — b'^) — b{b + a — 1)^ 

and show that the condition is valid for sufficiently great a. The domain in 
which condition (56) is not satished is bounded with respect to a. Since it 
is considered & < 4, a > 0, system (56) is equivalent to the following system 

b{b + a — 1)^ — 4a{b + ab — 6^) 

3(j2 

b{b + a — ly — Aa{b + ab — b‘^) 
a^(4 — b) 


{ r — 1 > 
r — 1 > 




r — 1 > 


r — 1 > 


dJ+i-h'-4(»+'’-?) 

3 

d^ + i-h'-4(»+'’-?) 


4-6 


Let be 


(57) 


= I - + 1 - - 


1,2 —2 + 6 ( 7 , — 6(7 + 1 — 3 ( 7 ^ 

= —P + - ^—b + -^- 


a 


a 


a^ a^ 

Then system (57) can be represented as 


(7^ 


(58) 

(59) 


r-l> -f{b,a) 


r — 1 > 


4-6 


f{b,a) 


(60) 


The function f{b,a) is quadratic polynomial in 6 with a positive coefficient 
of 6^. The abscissa of parabola peak is equal to 6o = 1 — 3a. If we consider 
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(T > |, then 6o < 0 and, therefore, on the interval b G (0,4) the maximum of 
function f{b, a) is attained for 6 = 4 and is equal to 




-3 + 


18 

— + 
a 


9 


a 


2 ■ 


For a > 3 + vT2 the value /(4, a) < 0. 

Since 6 G (0,4) and r > 1, for sufficiently great a we have 

J r - 1 > 0 > ^/(6,a) 

[r-1 >0> 


(61) 


i.e. condition (56) is satished. 

II. Let us show that for sufficiently great a and r > 1 the inequality 


ar > {b — a) (6 — 1) 


(62) 


is satished. Then the domain in which it is not satished is bounded on a. 
Condition (62) is equivalent to 

1 1 

r > —b'^-\ - - -6 + 1. (63) 

a a 

The right-hand side of the inequality is quadratic polynomial in 6 with a 
positive coefficient of 6^. The abscissa of parabola peak is equal to 6o = | + 

If we consider a > 7, then 6q > 4. In this case the maximum of the right- 
hand side of inequality for 6 G (0,4) is smaller than the value for 6 = 0, i.e. 
1. Thus, for sufficiently great a, 6 G (0,4) and r > 1 the relation 

1 o —a — 1, 

r > 1 > —b‘^ H- 6 + 1 

a a 

is satished, i.e. (62) is valid. 

III. If r > 1, 6 G (0,4), and a > 7, then condition (56) and inequality 
(62) are satished, i.e. the conditions of Theorem 3 are satished. □ 

Proof. (Remark 3) This is obvious from conditions (56). □ 
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